Symmetry-conserving vortex clusters in small rotating clouds of ultracold bosons 
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The properties of a special class of correlated many-body wave functions, named rotating vortex 
clusters (RVCs), that preserve the total angular momentum of a small cloud of trapped rotating 
bosons are investigated. They have lower energy and provide a superior description for the forma- 
tion of vortices compared to the mean-field Gross-Pitaevskii (GP) states that break the rotational 
symmetry. The GP vortex states are shown to be wave packets composed of such RVC states. 
Our results suggest that, for a small number of bosons, the physics is different from that of ideal 
Bose-Einstein condensates which characterize larger assemblies. 

PACS numbers: 03.75.Hh, 03.75.Lm, 03.75.Nt 



Mean-field descriptions of the many-body problem ex- 
hibit a ubiquitous symmetry-breaking behavior that ex- 
tends across several fields of physics, from nuclear physics 
[T], quantum chemistry [2], and metallic microclusters [3 
to semiconductor quantum dots j4| and trapped ultracold 
atoms HISIIS- 

Mean-field broken-symmetry solutions are expected to 
play the role of an effective ground state in the thermody- 
namic limit, N ~f oo, when quantum fluctuations about 
the mean- field state may be omitted |7j. A prominent 
example of such broken-symmetry states is given by the 
Gross-Pitaevskii (GP) vortex states in a harmonic trap 
[5] . Specifically, although each individual vortex carries a 
quantized amount of angular momentum [9J , the GP vor- 
tex solutions as a whole break the rotational symmetry 
of the confining harmonic trap |10j , and thus they are not 
eigenstates of the total angular momentum L = X]i=i '*■ 
In agreement with the general ideas of Ref. [7] , GP vortex 
states have been observed experimentally (see, e.g., Refs. 
[TTJ [T^J [131 mj ) for rotating Bose gases with large N. In- 
deed the energy advantage of symmetry restored states 
(see below) over the mean-field solutions diminishes as 
N increases (see section 1.2. in Ref. [4 ). 

Recently, the availability of optical lattices [T21 [T5] 
with a small number of particles per lattice site serves 
to motivate studies of small clouds of rotating bosons. 
For a small number N of atoms, however, quantum fiuc- 
tuations cannot be neglected, and one needs to consider 
methods beyond the mean-field approximation [3]. 

A natural way for accounting for quantum correlations 
about the GP vortex solutions in a harmonic trap is the 
method of restoration of rotational symmetry via projec- 
tion techniques. This method [4 was introduced recently 
in quantum dots [171 HH] and harmonic traps [3 [19] to 
describe individual particle localization and formation of 
rotating electron molecules (REMs) [20] and rotating bo- 
son molecules (RBMs), respectively. Here, we use pro- 
jection techniques to define and study vortex states with 
good total angular momentum, showing that these states 
can be properly referred to as rotating vortex clusters 
(RVCs). We stress again that the RVCs are eigenstates 



of the total angular momentum in contrast to the GP 
vortex states. For small N, the rotating-vortex-cluster 
states are the natural entities to be employed (in place 
of the GP vortices) for comparisons with exact solutions, 
which are eigenstates of the total angular momentum by 
their very nature. Furthermore we note the generality 
of the methodology of symmetry restoration. It applies 
as well to other broken symmetries, like spin symmetries 

We derive the RVC wave function by using an adap- 
tation of the two-step many-body method of symmetry 
breaking/symmetry restoration. We start with the ob- 
servation that the many-body GP vortex solution, ^'~^^, 
(as well as any mean- field solution exhibiting a breaking 
of the rotational symmetry) is a wave packet, and thus 
it can be expanded as a linear superposition over eigen- 
states ^N,L of the many-body Hamiltonian TC with good 
total angulal mommentum L, i.e.. 
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(ri,r2,...,rAr) = ^ CL$jv,L(ri,r2, . . . ,rAf). (1) 



For the two-dimensional case considered here (which 
is appropriate for a rapidly rotating harmonic trap), the 
eigenstates <&jv,l's can be approximated by using the pro- 
jection operator 
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which projects states with good total L out of the GP 
vortex state. The RVC wave functions are then given by 
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where |^]^ (6*)) is the original many-body GP vortex 
solution rotated by an azimuthal angle 9. The projec- 
tion techniques use the fact that the broken-rotational- 
symmetry states form a manifold of energy degenerate 
states (i.e., their total energy is independent of the az- 
imuthal angle 9); in this respect, the phases e'^^ in Eq. 



(l3| are the characters of the rotational group in two di- 
mensions. 

The expansion coefficients C^ in Eq. (fTl, which 
specify the spectral decomposition of the many-body 
Gross-Pit aevskii vortex state, can be calculated us- 
ing the projected wave functions (pi) for the ^a^.l's in 
the r.h.s of Eq. (fTl). Taking into consideration that 

^^n) = Hill </'o(ri)> one finds 
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where the overlap kernel is given by n{9) = 
{(j)o{0 = 0)\(j)o{S))^ , and the multiply occupied single or- 
bital (/>o(r) is a solution of the familiar Gross-Pitaevskii 
equation 

[H{v) + g{N - l)|(/.o(r)|2]c/.o(r) = eo0o(r), (5) 

with the single-particle Hamiltonian given by H(r) = 
p^/(2m) ~ 0.1 + muj^Y^ /2, where Q, is the rotational fre- 
quency of the trap and ljq characterizes the circular har- 
monic confinement. 

The total energy of the RVC is given by 



pRVC 



h{e)e"''^de 



n{0)e""^de, (6) 




with the Hamiltonian kernel being h(6) 



0)|7i|^}^ (6*)), where the many-body Hamil- 



tonian is Ti 
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above constitutes an effective continuous configuration- 
interaction scheme which lowers the mean-field energy by 
introducing correlations [? . The lowering of the ground- 
state energy brought about by the angular-momentum 
projection can be seen (see Ref. [21]) from evaluation 
of the GP ground-state energy using the spectral de- 
composition given in Eq. M. This yields the expression 
S^^ = El ICLl'Ef;^ with El \Cl\' = 1. Since E^^ 
is expressed as a weighted average of E^[' with positive 
weights, it is obvious that at least one of these energies 
obeys Ef^^ < E^^. 

To illustrate the essential qualitative difference be- 
tween the RVCs and the GP vortex states, we contrast 
in Fig. [I] their single-particle densitites (SPDs) for the 
case of TV = 9 trapped bosons and when the GP solu- 
tions exhibit either a (0,4) single-polygonal-ring of four 
[Fig. flja)] or a (2, 5) double-polygonal-ring configuration 
of seven [Fig. [ijb)] locahzed vortices (for the trap and 
other parameters employed, see the caption of Fig. [l]); 
(ni, n2) denotes rii (712) vortices on the inner (outer) ring. 
In sharp contrast to the GP single-particle densities, the 
RVC SPDs are circularly symmetric: the one [Fig. flic)] 
corresponding (through the aforementioned projection) 
to the four GP vortices exhibits instead a single contin- 
uous ring of depleted matter, while the other one [Fig. 



FIG. 1: Rotating vortex cluster and GP vortex solutions for 
N = 9 trapped bosons in a rotating trap with angular frequen- 
cies of trap rotation Sl/o^o = 0.5 [top three panels in the left 
column] and Q/uq — 0.565 [all four panels in the right column 
plus (h)] . For Q/uio = 0.5, a (0, 4) single polygonal ring of four 
vortices is involved, while for Q,/ujo — 0.565 a (2, 5) double 
polygonal ring of seven vortices develops. (a,b) GP single- 
particle densities (SPDs). (c,d) RVC single-particle densi- 
ties, (e-h) RVC conditional probability distributions (CPDs), 
with the fixed point (marked by a thick vertical arrow) at 
ro = (0,2.07Zo) (e), ro = (0, 1.14io) (f), and ro = (0, 3.2Zo) 
(g,h). (h) A horizontal slice of the CPD in (g) which mag- 
nifies the five vortices of the outer ring. The RVC total an- 
gular momenta are L = 28 [(c) and (e)] and L — 36 [(d), (f- 
h)]. The corresponding GP total-angular-momentum aver- 
ages ({*^^|L|*^'')) are 26.93 and 38.11, respectively. Note 
the elliptic shape of the vortex cores in the CPDs [see (e- 
h)] reflecting azimuthal fluctuations in the RVC state. The 
strength of the interparticle repulsion was taken Rg — 50 
[Rs = gm/ {2-Kfi^), see Ref [5j]. Unit length: Iq = ^ h / {rnujo) . 



md)] corresponding to the seven GP vortices exhibits in- 
stead two concentric continuous rings of depleted matter. 
Due to the symmetry restoration, the vortex struc- 
tures become "hidden" in the RVC single-particle den- 
sities. However, they can be revealed through the use of 
conditional probability distributions (CPDs) defined as 



P(r,ro) = ($ 
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FIG. 2; Spectral decomposition [Cl coeffiecients modulus 
square, Eq. Q] for GP vortex ground states in characteristic 
cases. Angular frequency of the trap; (a) n/uo ~ 0. (b) 
fl/ujo = 0.2. (c) n/tjo = 0.45. (d) n/coo = 0.55. The total 
angular momenta associated with the largest coefficients are 
Lmax ~ 0, 12, 21, and 36, respectively. The polygonal ring 
configurations of the GP vortices are also marked as (721,712). 
The ratio of the interparticle repulsion and the kinetic energy 
was taken Rs = 50 (see caption of Fig. [T|. In the figure, 
n is given in units of ujq- The GP total-angular- momentum 
expectation values are (a) 0.0, (b) 11.01, (c) 21.42, and (d) 
36.53. 



CPDs give the probability of finding a boson at position 
r given that another boson is located at a fixed point Tq. 

In the CPDs calculated for the RVC states [see Figs. 
flTe,f,g,h)], the fixed point is associated with a hump (lo- 
cal niaxiniuni) and the vortices are given by depressions 
(local minima) in the matter density. The number of 
vortices of a GP state and of an RVC projected out from 
it is the same for all fi's; e.g., in the CPD in Fig. W[c) 
four vortices are seen corresponding to the four GP vor- 
tices in Fig. [IT a). This refiects the fact that the projec- 
tion maintains the same (intrinsic or hidden) point-group 
symmetry. 

For the (2, 5) double-ring RVC, the fixed point can be 
placed on the inner or the outer ring. In the first case, the 
calculated CPD [Fig. flTf)] shows two vortices on the in- 
ner ring and remains uniform along the outer ring, while 
in the second case the calculated CPD [Fig.[IJg,h)] shows 
five vortices on the outer ring and remains uniform along 
the inner ring. This suggests that the rings rotate inde- 
pendently of each other in analogy with the case of rotat- 
ing boson molecules [12] and rotating electron molecules 
in quantum dots @]. 

To further investigate the wave-packet properties of the 
GP vortices, we have numerically calculated their spec- 
tral decomposition in terms of RVC states [see Eq. (fll]. 
In Fig. [2I we plot the expansion coefficients C^ [calcu- 
lated numerically according to Eq. (H] for A^ = 9 bosons 
and for several characteristic angular frequencies. The 



n = case is a rather trivial one where the GP solution 
preserves the circular symmetry, exhibits no vortices, and 
coincides with the corresponding RVC for L = (in this 
limiting case, one has Cq = 1 and C^ = for any L ^ 0). 

For the chosen value Rs — 50 (see the caption of 
Fig. fTl, non-trivial cases arise for fl/ujQ > 0.175. For 
n/wo = 0.2 [Fig. ^h)] and n/too = 0.45 [Fig. ^c)] 
the GP vortex solutions exhibit a (0,2) and (0,3) sin- 
gle polygonal-ring configurations, respectively. The ex- 
pansion coefficients Cl clearly demonstrate that the GP 
vortex states given as examples in these figures can be re- 
constructed from linear superpositions of RVC states [ob- 
serve the many non-vanishing values of Cl in Fig. 0^b) 
and Fig. [2Fc)]. Of central importance is the selection 
rules that the RVC angular momenta must obey in order 
to be compatible with the (ni,n2) intrinsic RVC point- 
group symmetry (which coincides with the explicit point- 
group symmetry of the associated GP vortex states) . In- 
deed, for the (0,2) vortex ring the RVC angular momenta 
obey the relation L = 2k, while for the (0,3) case, one has 
L = 3k, with fc = 0, 1, 2, .... The RVC angular momenta 
vary in a stepwise manner, and the value of the step coin- 
cides with the number of GP vortices on the single ring. 
For il/ujQ = 0.55, the GP vortex state exhibits a double 
(2,5) polygonal-ring structure, which imposes upon the 
RVC decomposition a more complex, but approximate, 
selection rule L = 2fci -I- 5fc2, with both ki and fc2 be- 
ing positive integer numbers. While the larger Cl's [Fig. 
l2jd)] conform to this rule, there are several smaller Cl's 
whose angular momenta fall outside this rule. The reason 
is that in the GP solution [see Fig.llTb)] the arrangement 
of vortices on the outer ring deviates slightly from being 
a perfect regular pentagon [22 . 

From the many RVC states that contribute at a given 
il to the spectral decomposition of the GP vortex states 
(including in the latter both ground and low-lying excited 
configurations), there is one with lowest energy Eqo^'^, 
which is the ground state within the RVC approxima- 
tion at this specific rotational frequency. For the fi's 
and the parameters considered in Fig. [2] we found that 
the RVC ground states have angular momenta Lqs as- 
sociated with the largest coefficients \Cl\ in the GP de- 
compositions; these Lgs's always obey the polygonal-ring 
selection rules discussed above. 

As aforementioned from the general theory of projec- 
tion techniques [2T], the RVC ground-state energies are 
lower than (or equal at most to) the corresponding GP 
ones for all values of the rotational frequency (Fig. l3|. 
Furthermore, the angular momenta associated with the 
RVC ground states are quantized and thus exhibit a step- 
wise increase as a function of the rotational frequency 
(see Fig. |4]). The magnitude of the steps in the ground- 
state RVC angular momenta changes with fi, since RVCs 
with different (n. 1,712) intrinsic vortex configurations be- 
come the ground-state as Cl increases. This behavior con- 
trasts with that of the ground-state GP angular momenta 
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FIG. 3: RVC (solid) and GP (dotted) ground-state energies 
for N = 9 bosons and Rs — 50, plotted versus the rotational 
frequency Sl/wo. The ground-state angular momenta are de- 
noted under the RVC curve. 



that vary continuously as a function of fi without any di- 
rect association to the (711,712) vortex configuration (Fig. 
111. In several instances, the RVC ground state has an 
intrinsic (711,712) point-group symmetry that is different 
from that of the GP vortex ground state at the same fi. 
This happens when the projection of an excited GP state 
results in a larger energy gain compared with that of the 
GP ground state. 

While we focus here on obtaining a proper symmetry- 
conserving vortex theory for a finite (small) number of 
trapped bosons, we comment on the RVC behavior com- 
pared with that obtained through exact diagonalization 
(EXD) calculations, which have become in recent years 
computationally feasible for smaller N. In particular, 
unlike the RVC case studied here (see Fig. |4]), EXD cal- 
culations in the lowest Landau level for A^ = 9 bosons ex- 
hibit quantized ground-state "magic" angular momenta 
that follow a Lj^ — niki + 7i2fc2 selection rule, but with 
the additional condition rii + n2 = N (unlike the RVC 
behavior where ni+n2 — q y^ N); see, e.g.. Fig. 10 in Ref. 
[23] and Fig. 2 in Ref. [24]. Moreover, it was shown in 
Ref. [25 that EXD ground-state wave functions describe 
formation of a rotating boson molecule [HI H] exhibiting 
two distinct (1,8) and (2,7) isomers of localized bosons. 

The above characteristic difference between the RVC 
and RBM [THJ US] solutions maintains for other values of 
N and Rs, reflecting the intrinsic point group symme- 
try of the angular momentum conserving theory (RVC, 
RBM via projection methods [l^ or EXD ^). This is 
particularly the case for low N and high Rs, where the 
RBMs are energetically favored and the RVCs may be re- 
garded as higher lying excited states. However, for higher 
N (occurring earlier for low Rs) the RVC may compete 
effectively with the localized RBM states, and eventually 
become the ground state. 

In conclusion, we have introduced a correlated many- 
body wave function, referred to as a rotating vortex clus- 
ter, which conserves the total angular momentum and 
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FIG. 4: RVC (solid) and GP (dotted) ground-state angu- 
lar momenta for A'' = 9 bosons and Rs = 50, plotted versus 
the rotational frequency Q/uJo- The quantized RVC angu- 
lar momenta are explicitly denoted in the figure, along with 
the corresponding polygonal-ring configurations (rii, 712). The 
Gross-Pitaevskii (ni, 712) configurations are also denoted. The 
GP angular momenta are not quantized; they are given by the 
expectation values {*^^|L|*5v^). 



has lower energy compared to the Gross-Pitaevskii solu- 
tion. The RVC is better suited to describe formation of 
vortices in small rotating clouds of trapped bosons com- 
pared to the mean-field GP vortex states that break the 
rotational symmetry. The GP vortex states were shown 
to be wave packets composed of such RVC states. The 
calculation of the properties of rotating-vortex-cluster 
states allowed for comparisons of qualitative signatures 
(e.g., ground-state angular momenta sequences) between 
the RVC and exact-diagonalization results. We con- 
clude that the physics of small rotating bosonic clouds is 
markedly different from that of larger assemblies known 
to behave as ideal Bose-Einstein condensates (properly 
described by the broken symmetry GP vortex solutions) . 
We hope that these results will motivate further experi- 
mental research in the area of correlated states in small 
bosonic systems (see, e.g., Ref. [IS]). 

This work was supported by the US D.O.E. (Grant No. 
FG05-86ER45234). 
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